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Abstract 

Spontaneous phase separation instabilities with the formation of various types of charge and 
spin pairing (pseudo)gaps in U > Hubbard model including the next nearest neighbor coupling 
are calculated with the emphasis on the two-dimensional (square) lattices generated by 8- and 
10-site Betts unit cells. The exact theory yields insights into the nature of quantum critical 
points, continuous transitions, dramatic phase separation instabilities and electron condensation 
in spatially inhomogeneous systems. The picture of coupled anti-parallel (singlet) spins and paired 
charged holes suggests full Bose condensation and coherent pairing in real space at zero temperature 
of electrons complied with the Bose-Einstein statistics. Separate pairing of charge and spin degrees 
at distinct condensation temperatures offers a new route to superconductivity different from the 
BCS scenario. The conditions for spin liquid behavior coexisting with unsaturated and saturated 
Nagaoka ferromagnetism due to spin-charge separation are established. The phase separation 
critical points and classical criticality found at zero and finite temperatures resemble a number 
of inhomogeneous, coherent and incoherent nanoscale phases seen near optimally doped high-T c 
cuprates, pnictides and CMR nanomaterials. 



PACS numbers: 65.80.+n, 73.22.-f, 71.27.+a, 71.30.+h, 75.10.Lp, 75.10.Jm 



I. INTRODUCTION 



A key element for understanding the complexity and perplexity in high-T c cuprates, 
manganites and colossal magnetoresistance (CMR) nanomaterials is the experimental ob- 
servation of phase separation (PS) instabilities at the nanoscale signaled by spin-charge 
separation and quantum phase transitions (QPTs) iMa). A new guiding principle for the 
search of new materials with enhanced Tc is the proximity to quantum critical points (QCPs 
for spontaneous first order QPTs attributed to intrinsic spatial inhomogeneities (see Ref. |6| 
and references therein). Strong quantum fluctuations dominate thermal fluctuations and 
affect the classical properties well above absolute zero temperature [7J. The inhomogeneous 
concentrated system in equilibrium can be well approximated as a quantum gas of decou- 
pled clusters, which do not interact directly but through the grand canonical ensemble, with 
different electron number per cluster. Our results for possible spatial inhomogeneities are 
directly applicable to nanoparticles and respective bulk nanomaterials which usually contain 
an immense number of isolated clusters in contact with a thermal reservoir by allowing elec- 



tron number per cluster to fluctuate. The finite-size optimized clusters may be one of t 
few solid grounds available to handle this challenging problem in a bottom-up approach 



lie 



;ap order parameters in the absence of a 
}]. The PS instabilities and spin-charge 



by defining canonical and grand canonical local 
long-range order, spin or charge density waves 
separation effects in bipartite Hubbard clusters driven by on-site Coulomb interaction U 
display QCPs which strongly depend on cluster topology 1CH13|. In frustrated (nonbipar- 
tite) geometries spontaneous transitions depend on the sign of the coupling t and can occur 
for all U by avoiding QCPs (level crossings) at finite U. The existence of the intrinsic QCPs 
and inhomogeneities associated with the PS instabilities, are crucial ingredients of the su- 
perconducting (SC) and ferromagnetic QPTs, providing important clues for understanding 
the incipient microscopic mechanisms of pairing instabilities in real space due to coexisting 
high/low electron (hole) or high/low spin up (down) densities in high-T c superconductors 
(HTSCs) and colossal magnetoresistive (CMR) nanomaterials respectively. However, small 
systems suffer from finite-size (edge) effects, so it is unclear whether the observed insta- 
bilities can survive in the thermodynamic limit. Thus, tests on reduced boundary effects 
are necessary to confirm the picture of local instabilities in larger systems in the so-called 
"optimized" Betts building blocks (finite square lattices) [ijj. A square infinite lattice is 
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FIG. 1: The 8-site finite unit cell (block) repeating periodically fill the entire (infinite) square 2d 
lattice. The cells can be defined by edge vectors l± = (2, 2) and h = (2, —2) (arrows in the figure) 
as defined in Ref. 11511. 



tiled by identical square unit cells containing L sites which periodically repeat in the lattice. 
For example, 8-site unit Betts' cell in Fig. [T]is used to fill the whole lattice. The square units 
restore basic symmetrical properties of the infinite square lattice and periodicity of clusters 
partially preserves translational and rotational symmetries of the original lattice. Therefore, 
Betts cells are believed to be the most preferred unit blocks for relieving frustrations over 
other structures with the same size. All the Betts unit cells can be defined uniquely by two 
edge vectors [l5| which represent translational vectors between the nearest two clusters. The 
lattice generated by Betts unit cells provides useful insights into certain physical aspects of 



the phase diagram in the t — J model 



161 ] . To our knowledge, an exact calculation of phase 



separation and pairing under doping has not been attempted in the Betts lattices applied to 
the Hubbard model either with nearest or next nearest neighbors. Different two-dimensional 
(2d) square structures, defined by the condition m 2 + n 2 = y/~L + \/~L with a linear size L (m, 
n are integers) , can be used as plaquettes to extrapolate the results to the infinite square lat- 
tice. Here our primary goal is an exact study of critical instabilities in the two-dimensional 
8- and 10-site Betts (generated) lattices. 



II. PHASE SEPARATION INSTABILITIES 



The lattice in two dimensions can be tiled by periodically repeated Betts (isotropic) unit 
cells with reduced edge (boundary) effects. 

In Fig. [T] an infinite square lattice has been tiled by L-site squares with edge vectors 
li = (2, 2) and I2 = (2, —2), which represent displacements of one vertex to the equivalent 
vertex in the neighboring tile with which it shares an edge. Notice that, if the sites are 
numbered as in Fig. [fl each odd site in the plaquette is surrounded by all the even sites (as 
nearest neighbors) and vice versa. The Betts unit cells take full advantage of the local space 
group symmetries of the isotropic 2d (square) bipartite lattice. The two-orbital Hubbard 
model with p bonding orbitals after elimination of the Cu sites can be reduced to a single- 
orbital Hubbard model with the nearest and next-nearest-neighbor (NNN) couplings by 
allowing holes to move within a given oxygen sublattice. Thus, we consider the minimal 
Hubbard model 

H = - l A Cja + U Y n n h lh 

where summation goes through all lattice sites L with coupling integral tij equal to t for 
the nearest and t n nn for the next nearest neighbors. The lattice frustration with t nnn 7^ 
allows study of electron pairing in the absence of electron-hole symmetry. Below, an exact 
diagonalization technique is used to extract the pairing instabilities and QCPs in finite 
8- and 10-site cluster-based Betts square lattices with periodic boundary conditions. In 
our previous work (Refs. |llMl3|). we have discussed PS instabilities in selected cluster 
geometries. A collection of "clusters" can be treated at a fixed average number of electrons 
(N) in a canonical ensemble or a fixed chemical potential /1 in a grand canonical ensemble. 
In the canonical ensemble we define a charge gap A C (N, T) at a given U and temperature T 
as A c (iV, T) — n + — n-, where //+ = E(N + 1, T) — E(N, T) and /i = E(N, T) - E(N - 1, T) 
energies are the first derivatives of the many-body (average) energy E(N,T) for adding or 
subtracting one electron in the charge sector for the iV-electron state. The canonical charge 
gap also describes electron fluctuations between different many-body cluster configurations, 
d +d N — > d N+1 + d N ~ 1 with different electron densities, n = N/L. The charge fluctuations 
closely resemble resonant valence fluctuations [l7| . Under certain conditions the charge gap 
can be associated with either an excitation (A c > 0) or a binding (A c < 0) energy. 



Physically, A c > manifests a stable d N state, while A c < implies a discontinu- 
ous spontaneous PS instability, i.e., first order transition (a negative compressibility) with 
the generation of d N+l and d N ~ x states of an inhomogeneous (mixture) of hole-rich and 
electron-rich phases coexisting with different n. Here, we focus on one hole off half fill- 
ing, (N) = L — 1. At the same time condition A c < implies derealization of holes 
(electrons) and pairing, which strongly suggests a possible local (correlation) coexistence of 
pairing with phase separation (segregation) rather than their competition. In the negative 
charge gap region (A c < 0), the spin pseudogap is nonnegative, i.e., A h > 0, since in the 
grand canonical ensemble we define it as the separation between the two consecutive peak 
positions of magnetic susceptibility Xs = dS z (/i,h)/dh for a given \i [9|]. However, when 
A c > (U > U c ) we define the canonical spin gap A s as an energy difference between 
the ground state with spin S and the lowest excited state with spin S' in the spin sector, 
A s = E(N, S',T) — E(N, S,T). When the excited state has a higher (lower) total spin, 
the spin gap is positive (negative). In addition, we introduced the grand canonical charge 
pseudogap A M > as apeak-to-peak distance between two consecutive peak positions of 
charge susceptibility jsHJ Xc = dN(fi)/dfi. A key element for the understanding of various 
electron instabilities is the exact relationship between the charge gap and pseudogap A c , 
A 1 ', and their corresponding spin counterparts A h , A s . The vanishing gaps and pseudogaps 
(nodes) are also crucial. As temperature approaches zero (T — > 0), the possible sign change 
in canonical gaps signifies the existence of QCPs related to first order (dramatic) changes, 
while zeros in the grand canonical pseudogaps, A M , A h , describe QCPs for continuous QPTs 
into quantum liquids. At finite temperatures, the vanishing gaps and pseudogaps define the 
classical critical temperatures (CCT) for pair condensations T c p , Tf and crossovers T c , T* 
into classical charge and spin liquids respectively. 

These pairing scenarios combined together give rise to the fascinating physics of Bose- 
Einstein condensation (BEC) and ferromagnetism due to competing high/low charge and 
high/low spin densities. The bound electron (hole) pairs with coupled anti-parallel (singlet) 
spins are complied with the Bose-Einstein statistics. Fig. [2] shows variations of gap order 
parameters A c , A M , A h , and A s with U in the 8-site Betts lattice at T — > 0, where we can 
see the number of QCPs and quantum phases. The broken line at < U < U c for a metallic 
gapless charge excitation (A^ — > 0) implies a quantum 'cold' charge liquid behavior. A finite 
non-zero field is required to break the coupled anti-parallel spins, A h > 0, which provides 
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FIG. 2: The charge A c and spin A s excitation gaps in the 8-site Betts lattice for (N) = 7 as 
a function of U at t nnn = and T — > 0. The charge and spin gap nodes, A C ' S (U) = 0, define 
a number of QCPs, U c , Ujjf, and Uf for various PS instabilities and complex phases. Phase I, 
< U < U c : a ground state describes coherent pairing (A h = — A c ) of electrons (holes) with 
anti-parallel spins and gapless charge excitations, A^ — > 0. Phase II for all U > U c : a stable charge 
ordered MH insulator, A M = A c > 0, coexists with a spin liquid, A h — > 0. At U c < U < Ujjf'- a 
spin liquid phase with gapless, low spin-3/2 excitations. At Uuf < U < Uf: a spin liquid phase 
with gapless, high spin-5/2 excitations coexisting with unsaturated ferromagnet. At U > Uf: a 
fully saturated spin-7/2 ferromagnet. Inset shows an enlarged area in the vicinity of key QCPs 
found for elementary square geometry in Ref. 

rigidity and stability for bound hole pairs (A c < 0). The equal (pseudo)gap amplitudes of 
paired charge liquid and coupled anti-parallel spins, A h = \A C \, called coherent pairing in 
Ref. pjj, are somewhat similar to conventional BCS-type pairing with a unique quasiparticle 
gap. 

This picture of coherent pairing of charge and spin entities in real space in Phase I (Fig. [2]) 
suggests a full BEC and possible superconductivity of electrons. In contrast, holes in Phase 
II, localized on separate clusters due to electron-hole pairing, A c = A M > 0, at U > U c behave 
as a low density Mott-Hubbard (MH) excitonic insulator. The 'cold' charge localization at 
the QCP, U c , gives rise to a crystallized nonconducting (Mott) state driven by correlation 
U (decreasing t or pressure, p). Hence, U c signifies a QCP for a continuous metal-insulator 




FIG. 3: The charge gap A c as function of U in the 10-site Betts lattice for (N) = 9 at t nnn = 
and T — > 0. In contrast to the 8-site lattice, optimized pairing is shifted to larger U values. At 
large U, there is a phase transition from a MH insulator into an unsaturated ferromagnetic ground 
state with QCP, Ujjf = 28.6. Inset displays the details in a moderate U region, where equal 
(pseudo)gaps, A h = — A c , display coherent pairing. 



(MI) transition from a charge liquid, A M — > 0, into a MH insulator. Notice that this QCP 
signals also spin-charge separation at U > U c [8]. The evolution of spin gaps A s and A h 
is also shown in Fig. [2j The isolated QCP, U c , signifies a 'cold' quantum melting of anti- 
parallel spin- insulator, A h > 0, into a degenerate spin-liquid without rigidity (A h — > 0). The 
quantum derealization of spin as U increases reveals an opposite trend to charge behavior. 
The nodes of the spin gap, denned by A S (U) = 0, are the QCPs for the onset of spin PS 
and spin density (nanophase) inhomogeneities due to spontaneous redistribution of the spin 
liquid among the clusters in the negative gap region. Magnetic field fluctuations lead to 
segregation of clusters into regions rich in spin ~f and spin J,, i.e., formation of magnetic 
domains. The QCPs for spontaneously broken symmetry at Ujjf = 14.9 and Up = 19.9 in 
Fig. [2] correspond to phase transitions into an unsaturated ferromagnet (Up) and a saturated 
Nagaoka-like ferromagnet [llj. Phase at U c < U < U\jf in Table [J is a degenerate spin-3/2 
liquid without spin rigidity (A h — > 0). The sign change in A s at Ujjf < U < Up implies the 
presence of a QCP between Uuf and Uf- The PS region A s < corresponds to unsaturated 
ferromagnetism, i.e., S = 5/2, while A s > region is a gapless S = 5/2 spin liquid (A h — > 0) 
with a finite spin excitation gap (A s > 0) for S = 7/2. 



TABLE I: The energy levels in the 8-site cluster for various S at N = 7. At Uc < U < Uuf, a 
spin-3/2 liquid (A h — > 0) has a spin-5/2 excitation gap, A s > 0. In Uuf < U < Uf range, the 
state with A s < is an unsaturated 5 = 5/2 ferromagnet, while next region with A s > for 
increasing 5 = 7/2 is a gapless spin-5/2 liquid. At U > Uf, a saturated ferromagnet (A s < 0) 
behaves as a spin-7/2 liquid. 
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FIG. 4: A contour plot of the canonical charge gap as a function of U and T in the 8-site Betts lattice 
at average (N) = 7. The negative gap region shrinks and totally disappears at about T > 0.16. 
The solid contour denotes a second order phase transition for classical critical temperature, (U), 
at which the gap vanishes, A C (U, T^) = 0. Phase I with a shaded area is the metallic charge pairing 
phase for spontaneous first order PS instabilities and density inhomogeneities. Phase II without 
a shaded pattern is the stable MH insulator region. The broken line for onset of coherent pair 
condensation of anti-parallel spins at Tf (U) is indicative of BEC and superconductivity. Below 
Tg the paired anti-parallel spin stiffness of bound charge describes coherent pairing phase, while 
above Tf preformed pairs with unpaired spins is characteristics of incoherent phase [111. Il9j|. 
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FIG. 5: The density (N) versus m in the 8-site Betts lattice at three U values near U c in the grand 
canonical ensemble at T — >• 0. The MH plateau near (N) = 7 at U = 10 vanishes with abrupt 
change at U = 4 and 8 QCPs, jXpiJJ). In the canonical ensemble at (N) = 7 the change from 
N = 8 to N = 6 with unstable N = 7 clusters suggests a negative compressibility for the first order 
PS and inhomogeneities in a spinodal region [121 ] . 



Phase U > £7^: the region with A s < describes a spin liquid (A' 1 — > 0) of fully polarized 
Smax = 7/2 Nagaoka-like ferromagnet coexisting with the MH-like insulator, A c > 0. In 
Fig. |3] the exact gaps in the 10-site Betts lattice reproduce the key features of spontaneous 
PS instabilities. The negative charge gap region defines a first order PS, while the positive 
gap region is a signature for continuous, second order transitions. The oscillations in the 
Betts 10-site lattice at small and moderate U are due to reduced square symmetry similar 

lJl2|. d the 



to the anisotropic 2 Hubbard clusters 



2x2 Hubbard cluster both have higher rotational and reflection symmetries. Figs. [2] and [3] 
are consistent with our phase diagrams in 2 x 2 and 2x4 geometries lfj, llll]. Next we 
show that the effect of QCPs is also tangible over a wide range of temperatures. Contour 
isolines along which the charge gap attains a constant value, are illustrated in Fig. H] for 
tnnn = 0. The contour map defines the charge gap isolines as a function of T and U. As 
temperature increases charge gap changes sign and favors electron-hole pairing above the 
critical temperature T c p defined by A c (U,Tf) 

(second order) transition for the onset of charge pairing below T c p [2|, 



0. The vanishing gap describes a smooth 



12 



13]. Phase at 



T > T c p is the MH insulator (A c > 0), while at low temperature charge paired (A c < 0) 




FIG. 6: Charge and zero field magnetic susceptibility peaks \c and grand canonical Xsi denoted 
by T c and T* in the 8-site Betts lattice at U = 8. The zero-field spin susceptibility peak T* 
terminates at Tf, while charge peak T c (/j,) for closed pseudogap, — > 0, terminates at the QCP, 
fip. The spin peaks separate the low (S = 3/2) and high spin (S = 5/2) phases above and below 
T* respectively. With the onset of strong spin (singlet) pairing below T p spin collective excitations 
become gapped. The charge and spin peaks T c and T* emerge and reconcile for the Fermi liquid 
behavior at the close vicinity of T p . The increase of the distance between T* and T c , peaks above 
T p implies a spin-charge separation [9(. The charge pseudogap region due to the splitting of T c 
into two branches above T p manifests electron hole (pairing) order broken by temperature. As 
in [l|3], the formation of incoherent charge pairs at T c p < 0.06 results in a charge pseudogap area 
above Tf. 



phase is a metallic charge liquid, with A M = (Fig. |6]). 

The continuous correlation driven MI transition, broken by electron hole pairing at 
T > T c p , is also somewhat different from the conventional T — p phase boundary between 
the ordered (quantum) phase at low and disordered (normal) phase at high temperatures, 
described by an order parameter, which is non-zero in the ordered phase and zero in the 
disordered phase. The charge gap shows pairing below T c p and proximity to electron- hole 
pairing A c > Oand MH pseudogap behavior (A M > 0) above T c p . However, unlike in the BCS 
theory, as T increases A c differs from spin gap A h , consistent with the existence in Fig. H] 
of two energy (gap order parameters or) scales and two distinct consecutive condensation 



temperatures for paired charges, T^, and bound anti-parallel spins, Tjf [4J, The spin 
rigidity for bound anti-parallel spins (A h > 0) below Tf manifests the onset of coherent 
pairing and possible superconductivity ay , 111] . In Fig. H]the incoherent pairing of bound 
charge A C (T) 7^ without rigidity A h (T) = exists at Tf < T < T c p . The grand canonical 
result for the (dramatic) abrupt doping dependent behavior of (N) at < U < U c in Fig. 
provides strong evidence for the existence of a QCP at equilibrium in \i space associated 
with a spontaneous first order PS. The first order transitions at U — 4 and 8 at T — > 
display the QCP, fip, at optimal doping. The visible plateau near (N) = 7 in Fig. [5] for 
A c > at U = 10 describes the rigidity of the MH insulating charge gap 11]. 

The canonical negative gap at T = displays also a PS spinodal [/i + ,/x_] region, a+ 
jXp±A c /2 for coexisting inhomogeneous N?6 and 8 clusters with electron redistribution |2|. 



3, 



11| . The effect of quantum criticality can be felt also at low T without even ever reaching the 



ground state. Fig. [6] shows the susceptibility peak crossover T c (/i) and T*(/x) temperatures 
versus /x for gapless charge and spin excitations in the grand canonical ensemble. The 
peak-to-peak separation for the two splitting branches of T c beneath T* in the underdoped 
regime at p, > pp defines the so-called pseudogap above Tf . The distinct pseudogap phase, 
fundamentally different from pairing, is a precursor to the change of the charge gap sign at 
low temperature, T < T c p [18]. We find that the pseudogap predominantly originates from 
(N) = 7 ( pla teau) phase at U c < U < Uup, which is stabilized as T, fi or U increases or p 



decreases 



191 ] . In many real materials, the contribution of the t nnn term (frustration) can 



be important for breaking the particle-hole symmetry in hole/electron doped systems. 

Thus, it is necessary to consider the more realistic case, t nnn 7^ 0. The charge gap 
contour plots for different t nnn and U values at T = are shown in Fig. [7] by tuning U 
and the next nearest coupling, t nnn . The condition A C (U, t nnn ) = defines the boundary 
tn nn (U) between positive and negative gaps, which separates smooth transitions from first 
order phase transitions. The t nnn < has a detrimental effect, while t nnn > enhances 



pairing. The PS region shown in inset in Fig^ 



3 at t nnn = t manifests a coherent pairing 



instability for all U. As for the tetrahedron [13], the vanishing PS region at t nn n = —t, is 
signaled by stable d 7 MH configurations with S ma x = 7/2 ferromagnetism. 




FIG. 7: The contour plots of charge gap A c as a function of U and t nnn in the 8-site Betts lattice. 
As t nnn > increases, the crossover point shifts to a larger U value (the inset displays the charge 
gap at t nnn = t = 1 in the enlarged area of U). As t nn n approaches t or —t, PS instability region 
increases or shrinks respectively, by avoiding for t nnn = — 1 of level crossings (QCPs) at finite U. 



III. SUMMARY 



The Betts lattices, generated by finite clusters, provide strong support for PS instabilities 
found in generic 2x2 and 2x4 clusters 0, n]] . Many complex critical phenomena, apparent 
in approximate treatments of "concentrated" inhomogeneous systems without long-range 
order are naturally reproduced in the grand canonical and canonical ensembles of indepen- 
dent clusters, which do interact thermally. The exact calculations in optimized clusters are 
a powerful tool for unveiling hidden generic QCPs of phase separation, 'hot' and 'cold' lo- 
calization and melting of charge and spin entities. The theory yields strong evidence for the 
existence of a QCP at optimal doping, /ip, responsible for a first order PS transition from 
a SC metal into a MH insulator hidden beneath the SC dome in the HTSCs 2), [j, [3]. The 
negative canonical gaps are sufficient for spontaneous PS and accompanied inhomogeneous 
redistribution of the electron charge or spin at energy level crossings. Small variations in 
model parameters or intrinsic symmetries of the unit block can lead to dramatic qualita- 
tive changes in the vicinity of QCPs, which control the physics of phase separation over 
a significant portion of the phase diagram. The key intrinsic QCPs, extracted from exact 
calculations, decipher the mystery of quantum and classical critical behavior seen in transi- 



tion metal oxides, cuprates, pnictides and CMRs nanomaterials under variation of external 
parameters and provide a basic scenario for Nagaoka-type ferromagnetism and correlation of 



electron pairing with phase separation 



20M22| . Building blocks spontaneously organized into 

novel 



ordered structures suggest a bottom-up key paradigm to control QCPs and fabricate o 
assembled nanostructures with new superconducting and magnetic properties jsl, 23I . 



24]. 
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